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Abstract 

We consider a multidimensional Ito process Y = (lt)te[o,T] with 
some unknown drift coefficient process bt and volatility coefficient 
a{Xt,6) with covariate process X = (A^t)tg[o,r]) the function cr{x,9) 
being known up to G 0. For this model we consider a change 
point problem for the parameter 9 in the volatility component. The 
change is supposed to occur at some point t* £ {0,T). Given dis- 
crete time observations from the process {X,Y), we propose quasi- 
maximum likelihood estimation of the change point. We present the 
rate of convergence of the change point estimator and the limit there- 
oms of aymptotically mixed type. 
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1 Introduction 



The problem of change point has been considered initially in the framework 
of independent and identically distributed data by many authors, see e.g. 
Hinkley (1971), Csorgo and Horvath (1997), Inclan and Tiao (1994). Re- 
cently, it naturally moved to context of time series analysis, see for example, 
Kim et al. (2000), Lee et al. (2003), Chen et al. (2005) and the papers cited 
therein. 

In fact, change point problems have originally arisen in the context of 
quality control, but the problem of abrupt changes in general arises in many 
contexts like epidemiology, rhythm analysis in electrocardiograms, seismic 
signal processing, study of archeological sites and financial markets. In par- 
ticular, in the analysis of financial time series, the knowledge of the change 
in the volatility structure of the process under consideration is of a certain 
interest . 

In this paper we deal with a change-point problem for the volatility of a 
process solution to a stochastic differential equation, when observations are 
collected at discrete times. The instant of the change in volatility regime 
is identified retrospectively by maximum likelihood method on the approx- 
imated likelihood. For continuous time observations of diffusion processes 
Lee et al. (2006) considered the change point estimation problem for the 
drift. In the present work we only assume regularity conditions on the drift 
process. De Gregorio and lacus (2008) considered a least squares approach 
following the lines of Bai (1994, 1997) of a simplified model also under dis- 
crete sampling while Song and Lee (2009) considered a CUSUM approach. 
Finally it should be noted that the problems of the change-point of drift for 
ergodic diffusion processes have been treated by Kutoyants (1994, 2004), but 
the asymptotics and the sampling schemes are different from this paper. 

The paper is organized as follows. Section [2] introduces the model of 
observation, the regularity conditions and some notation. Section [3] studies 
consistency and the rate of convergence of estimator of the change while 
asymptotic distributions are considered in Section HI A mixture of certain 
Wiener functionals appears as the limit of the likelihood ratio random field, 
and it characterizes the limit distribution of the change-point estimator. 
Those sections assume that consistent estimators of the volatility parameters 
are available. Section [5] presents some practical considerations and a proposal 
to obtain first stage estimators of the volatility parameters which allow to 
obtain all asymptotic properties stated in the previous sections. Finally, 
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Section presents some numerical analysis to asses the performance of the 
estimators. Tables are collected at the end of the paper. 



2 Estimator for the change-point of the volatil- 
ity 

Consider a dimensional Ito process described by the stochastic differential 
equation 

dYt = btdt + a{Xt, 9)dWt, t E [0, T], (1) 

where Wt is an r-dimensional standard Wiener process, on a stochastic basis, 
bt and Xt are vector valued progressively measurable processes, and a{x, 6) 
is a matrix valued function. 

We assume that there is the time t* across which the diffusion coefficient 
changes from (t(x, ^o) to a{x,6i). The change point t* G (0, T) is unknown 
and we want to estimate t* based on the observations sampled from the 
path of (X, F). The coefficient a{x,6) is assumed to be known up to the 
parameter 9, while ht is completely unknown and unobservable, therefore 
possibly depending on 9 and t* . 

The sample consists of {Xt^,Yt^), i = 0, l,...,n, where ti = ih for h = 
hn = T/n. The parameter space of 6* is a bounded domain in M"'", do > 1, 
and the parameter ^ is a nuisance in estimation of t*. Denote by 9* the true 
value of 9i for i = 0, 1. 

Let = 1^1 ~^ol- We will consider the following two different situations. 

(A) 9q and 9^ are fixed and do not depend on n. 

(B) 9q and 91 depend on n, and as n ^ oo, 9^ ^ 9* e Q, 'dn and 
n-d"^ —>■ oo. 

In Case (A), is a constant "do independent of n. 

We shall formulate the problem more precisely. It will be assumed that 
the process Y generating the data is an Ito process realized on a stochas- 
tic basis B = F,P) with filtration F = {J^t)t£[o,T], and satisfies the 
stochastic integral equation 




Yo + /J bjs + /J aiXs, 9*o)dWs for t e [0, t*) 
Yt^ + hds + 9l)dW, for t e [t*, T]. 
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Here Wt is an r- dimensional F- Wiener process on B, and bt, Xt and o"(x, 9) 
satisfy the conditions below. Let A" be a closed set in (possibly X = W^'^) 
and denote the modulus of continuity of a function f : I ^ W^^ by 

wj{6J)= sup \f{s)-f{t)\. 

s,tei, \s-~t\<5 

For matrices A = (aij) and B = (bij) of the same size, we write A^"^ = 
A^A, A[B] = YjijO-iikj = Tt{A^B), and the Euclidean norm of A by \A\ = 
Set S{x,9) = cr(a;, 6*)*^^. \i{A) denotes the minimum eigenvalue 
of a symmetric matrix A. 

[H]j (i) a{x,t) is a measurable function defined on x [0,T] satisfying 

(a) inf(^.,0)g;v'xe Ai(5'(x,6')) > 0, 

(b) derivatives d^a (0 < i < j + [do/2]) exist and those functions are 
continuous on x 6, 

(c) there exists a locally bounded function L : X x X x Q ^ M_|_ such 
that 

\a{x, 6) - a{x\ e)\ < L{x, x' , e)\x - x'\" {x, x' eX^OeQ) 
for some constant a > 0. 

(ii) {Xt)t£[o,T] is a progressively measurable process taking values in X 
such that 

as n ^ oo. 

(iii) {bt)t£[o,T] is a progressively measurable process taking values in M'^ 
such that {bt — &o)te[o,T] is locally bounded. 

Remark 1. The term "locally bounded" in [H]j (i) (c) means, as usual, 
being bounded on every compact set. The case where the drift bf changes its 
structure at time t*, or any time in force, is included in our context because 
bt admits jumps. The case of time dependent a is included by making Xt 
have argument t. Needless to say, if we set X or a. part of X as F, then our 
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model can express a system with feedback, in particular, a diffusion process. 
By [H]j (ii), t Xf is continuous a.s. Also, [H]j (ii) imposes a restriction 
on the rate 'dn- For example, when a = 1, for a Brownian motion X, it 
suffices that n-i^^/logn —> oo, due to Levy property. The additional [do/2] 
time differentiability to j is used only in Step (iii) of the proof of Theorem 
[U Therefore, it is possible to replace the range of £ to "0 < £ < j" under a 
condition that ensures the the Hajek-Renyi type estimate just before going 
to Inequality (jlj) below. 

Write AiY = Yt^ - Ft._^ and let 

[nt/T] n 
i=l i=[nt/T]+l 

where 

G0) = \ogdetSiXt^_,,9) + h~'SiXt^_„9r'[{A,Yr% 

Suppose that there exists an estimator 6^ for each 6^, k = 0,1. Each 
estimator is based on (Xt., l^.)j=o,i,.--,n and so depends on n. To make our 
discussion complete, in case 91 are known, we define 9k just as 9k = 9^. This 
article proposes 

in = aTgmm^^[QJ,^^nit;9o,9l) 

for the estimation of t*. More precisely, t„ is any measurable function of 
(XtJj=o,i,...,n satisfying 

^n{in;9o,9i) = mm ^n{t;9o,9i). 

tG[0,T] 

3 Rate of convergence 

We introduce identifiability conditions in order to ensure consistent estima- 
tion. In Case (A) we assume 

[A] P[S{X,,;9*) S{X,,;9l)] = l; 

In Case (B) we assume 
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[B] E{Xt*,9*) is positive-definite a.s., where 

E{x,9) = (Ti{{dg,.oS)S~\dsu.jS)S~'){x,9)] ° ,9 = 



\ / h,i2=l 

Remark 2. Since 9*) is the Hessian matrix of the nonnegative function 
Q{x,9*,9) ■.= TT(^S{x,9*y^S{x,9) - Id^ -logdet (^S{x,9*)-^S{x,t 

of 9 at 9*, ^*) is nonnegative-definite. 

The following property will be necessary to validate our estimating pro- 
cedure. 

[C] \9k - 9l\ = Op{'dn) as n ^ cx) for A; = 0, 1. 

In case the parameters are known, 9k should read and then Condition [C] 
requires nothing. Section [5] presents an example of estimator for 9k which 
satisfies Condition [C]. 

Here we state the result on the rate of convergence of our change-point 
estimator. 

Theorem 1. The family — t*)}neN 'is tight under any one of the 

following conditions. 

(a) [H]i, [A] and [C] hold in Case (A). 

(b) [H]2, [B] and [C] hold in Case (B). 

In Case (B), this result gives consistency of t„ since nd^^ ^ oo by as- 
sumption. 

The rest of this section will be devoted to the proof of Theorem [T] Define 
a stopping time r = t{K) by 

t{K) = inf |t; \Xt\ + \bt\ > A'j A T 

for K > 0. X'^ denotes the process X stopped at r. Write Si{9) = S{Xl., 9), 
and AiY^ = Y.^ - Y.^ . Let 

[nt/T] n 
i=l i=[nt/T]+l 
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where 

gm = l{.>o}|logdet5,„i(^) + /i-iS,„i(^)-i[(A,n''']} 

= l{.>o}logdet5,„i(^) + /i-i5,_i(^^)"M(A.F-)®']- 
Then sup^g^ \9i{^)\ ^ for any compact set /C in 6 under [H]i. Denote 
by -Ejj^i the conditional expectation with respect to J^ti^^ under the true 
distribution for > t* . 

Lemma 1. Fort > t* , 

^n(t;^0,^l)-*n(t*;^0,^l) = M„(t;^o,^l)+^n(t;^0,^l)+Pn(t;^0,^l), 

where 

[nt/T] 

M^{t-9o,9,) = {[9^W-g^m-E[^M9o)-gM]}, 



An{t; 60, 61 



i=[nt*/T]+l 

it/T] 

E 

4nt*/T] 



j=[nf/T]+l ^ ^ ^ 

-logdet (^5,„i(eo)~'5.-i(^i))}, 



L'"7 ^ J f 

Pn{t;eM = l|.>o} 5^ Trj (5,_i(^i)-i-5,_i(^o)"') 

j=[nf*/T]+l ^ 

The proof of Lemma [1] is straight forward and omitted. 

Remark 3. Later we will consider substitution of estimators 6k to 6k, k = 
0, 1. Then the expectation E^^^[gi{6Q) — gi{6i)] is taken before the substitu- 
tion, and so 

[nt/T] . 

M4t; 60, k) = [gm - gi{6^)] - E^^{g,{6,) - g,{6^)] , , I 

i=lnt*/T]+i ^ eo=eo,e,=e,) 

In particular, the second term in the braces is not necessarily ^jj_i-measurable. 
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We will need a uniform Hajek-Renyi inequality. Let D be a bounded open 
set in W^. The Sobolev norm is denoted by 



s,p 



p 

LP{D) 



i=0 



1/p 



for / e W'^'P{D), the Sobolev space with indices {s,p). Suppose that p > 1 
and s > d/p. The embedding inequality is the following 



sup 1/(^)1 <C||/||,,, (feW^'^m 



9eD 



(2) 



where C is a constant depending only on s,p and D. We will apply this 
inequality for / G C^{D), and the validity of such an inequality depends on 
the regularity of the boundary of D; see e.g. Yoshida (2005) for the relation 
to the GRR inequality. 

Lemma 2. Let (fi, JF, F = {J^j)j^z+, P) be a stochastic basis. Let D be a 
bounded domain in M'^ admitting Sobolev's inequality ([2]) for some p G (1,2] 
and s G N such that s > d/p. Let (cj)jgz+ be a nondecreasing sequence of 
positive numbers. Let X = {Xj)j^i^ be a sequence of random fields on D for 
j G Z+ satisfying the following conditions: 

(i) For each {w,j) E^lx Z+, Xj G C'{D); 

(ii) For each {0,i) G -D x {0,1,. {dlXj{9))j^z+ is a zero-mean L^- 

martingale with respect to F. 

Then there exists a constant C depending only on s,p and D, not depending 
on X, such that 



1 

max — sup \Xj 

j<n Cj D 



for all a > and n G Z+. 



> a 



< 



Proof. Let B = Lp{D), then B is p-uniformly smooth; see Example 2.2 of 
Woyczyhski (1975), p. 247. We apply Theorem in Shixin (1997) to conclude 



max — \\dlXj\\B > 



j<n c 



< 



n 



j=0 J 



for i G {0, 1, s} for some constant Ci. Therefore yeilds the result. □ 

Proof of Theorem Ui For the proof, we may assume T = 1 for notational 
simphcity without loss of generality, 
(i) Let e be an arbitrary positive number. Set 



in Case (A), and set H{x) = \i(E{x,9*)) in Case (B). We denote (j(t;9) = 
o'{Xl, 6) and h{t) = H{Xl) in what follows. Those processes depend on K 
by definition while it is suppressed from the symbols. Set Bk = {r = 1} and 
fix a sufficiently large K so that -P[-Bx] < 

We notice that h{s) > and that h{t*) > a.s. on from the identi- 
fiability condition [A]/[B] since X^, = Xt* on Bx- We will show that there 
exists a positive constant q such that 



1-2 




< e. 



Define the event As by 




for 5 G (0, 1 - r) 



On Asi it holds that 




and also that, for t G [t* + 5, 1], 
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hit*). 



> 



2{l-t*) 



9 



Choose a 5 so that P[^5] > 1 — e/2 by the continuity of h, and next 
choose a positive number = c(e, 6) such that 
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2(1 -t* 



-h{t*) > 5q 



> P 



> 1 



6 



Then 
P 



1 f 

inf / h(s) ds < 5c. 



< e. 



As, inf - 

te[f,i] t 



1 /■* 

J h{s) ds < 5c, 



(ii) With Lemma [H we decompose \l/„(t; 9q, 9i) — \l/„(t*; 9q, 9i) as follows: 

■^r^{t■Jo,9,)-'^nit*■,9o,9l) = Mn{t;9o,9i)+Anit;9o,9i)+p^it;9oJi). 
Let M > 1. We have 



P[n^l{L -t*)> M] < P 



< P 



inf $n(t;^o,^i)<$n(r;^o,^i) 

t:n-al{t-t*)>M 

inf ^„(t;^"o,^i) <^n(t*;^o,^i) 

t:n-&'i{t-t*)>M 



where 



Pl,n 
P2,n 
P'i,n 



P 
P 
P 



sup 

inf 

\{.t-t*] 

sup 



t:n{j 



P2,n 


+ P3,n 




1 


[nt] 


- [nt*] 




1 


[nt] 


- [nt*] 




1 


[nt] 


- [nt*] 



P{Bk] 
(3) 



> 



Mn{t;9o,9i ^ 

An{t;9o,9,)<cX 
Pn{t;9o,9i) > 



c d'^ 



C 'ff^ 



Here we read inf = oo and sup0 = — oo. We will estimate these terms, 
(iii) Estimate of Pi^„. In Case (B), let 



[nt] 

Mn{t;9) = {^eg^{9) - EZ^[^eg^{9)]] 

i=[nt*]+l 
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Let be an open ball such that 6* E Q and C G. Since 



sup |M„(t; ^0,^0 1 1^0-^1 r' < snpMnit; 



one has 



sup 



.t:n'&lit-t')>M [nt] - [nt*] 
+P[\0o -0i\> 2?9„] + P[^o ^ e] + P[^i ^ 6] 



6 



< P 



sup 



t:ni)l{t-t*)>Ai [nt] - [nt*] 



snp \Mn{t;9)\ > 



6 



+P[\9o -9^\> 2i9„] + P[9o ^ 6] + P[9^ ^ 6]. 

By the uniform version of the Hajek-Renyi inequality in Lemma [2] applied 
to the case p = 2, s = 2 + [(io/2] and D = Q, we see under [H]2 that 



P 



sup 



sup \Mn{t;9)\ > 



< 



C 



cm 



t:n^l{t^f)>M [nt] - [nt*] ' ' " 5 

therefore 

lim„_»oo-Pl,n < Pe(M) (4) 

thanks to 

P[|^o-^i| >2^„] < P[\9,-9l\>U^]+P[\9,-9l\>U^] 
for large n. 

In Case (A), Let 6^ be an open ball such that 9^ C 6 and G Ofc for 
each /c = 0, 1. 



1 



sup 



t:n'dl{t-t*)>M [nt] — [nt*] eggeo 



sup |M„(t;^o,^i)| > 



+P[9o^Q^]+P[9i^Qi]. 
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We apply the Hajek-Renyi inequality for Mn(t;6Q,6i), which is a difference 
of two random fields on 6^ to be done with one by one, in order to obtain 
dH) under [H]i. 

(iv) Estimation of P2,n- First we consider Case (B). There is a positive 
constant C2 independent of n such that 

Tr(^5,„i(^o)~'5.-i(^i) -/d) -logdet (^5,_i(^o)"'5.-i(^i)) 

> {Ai(S(X,^_^,r))+r„,_i}|^i-^;p 
for all i, where maxj iT-n^j-il < C2'&n, on the event 

BK,n = BKn{eo,eieQ,\h-e*\<^^{k = o,i)}. 

Thus 

1 



P2,n < P 



inf 



t:n-dl(t-t*)>M [nt] — [nt*] 



^0,^1) 1^1 -^or'<4Q, Bk, 



-p 



< p 



inf 



+ P[Bkn] 



[nt] 



t:n-&l{t-f)>M [nt] — [nt*] 



J2 {Xi{E{Xl_^,e*))+r^,^,}<Ac, 



--[nfj+l 



for large n. The scaled summation converges to the corresponding scaled 
integral uniformly in t a.s., hence from Step (i) we have 



lim„^ooA,n < P 



inf 



1 



te[t*,i] t -t* 



h{s) ds < 5ce 



< 2e 



for large n. 
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We will consider Case (A). There is a positive constant C2 independent 
of n such that 



TrLS,_i(^o)~'5,„i(^i) -lA- logdet 5,_i(^^o)"'5^-i(^i) 



> ii\^s,^^{eir^Si^i{ei) - hj - logdet (Si.i{e*y's,.i{ei) 
-C2i\0i-e;\ + \eo-e;\) 

for all i on the event -B^„ = Bk H {^0 £ ©O) ^1 G 9i} because there exists a 
continuous derivative dga by [H]i. In this way, 



P2,n < P 



inf 



t:n(t-t*)>M [nt] — [nt*] 



Therefore, 



lim„^oo^2,n < P 

< 2e 

by Step (i). 

(v) Estimation of P^^n- We have 



inf 

te[t*,i] t-t* 



h{s) ds < 5ce 



+ e 



sup 



S{Xt,9k)-S{Xt,9l 



<C^n (A: = 0,1), 



sup 

t€[t*,l] 

and 



S{XtJ,r'-S{Xt,9, 



k) 



1 



{|6lfe-9*|<2,?„}nBK 



<C'&n (A; = 0,1) 



sup 

i:>[nt*]+2 



n 



In the last estimate, the local a-Holder continuity of a was used. Then on 

BK^{\ek-ei\ <2^„ (fc = o,i)}. 



sup - — — - 

t:ni)l{t-t*)>M [nt\ - [nt*\ 



(5) 
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because of [H]j (ii). Consequently, we see lim„^oo-P3,n < e due to [C] and the 
localization by B^- 

(vi) From the estimates in Steps (ii)-(iv) and making K sufficiently large, we 
have 

]hRn^^P[ni^l{L - t*) > M] < p,(M)+5e 
for any M > 1 and e > 0. Therefore, 
limM^oolim n—>oo 

which shows the tightness of — In a quite similar way, we can 

show that {ri'dl^itn — t*)-}n is tight, and hence the family {n{}'^(tn — t*)}n is 
tight. □ 

4 Asymptotic distribution of the change point 
estimator 

This section discusses limit theorems for the distributions of the estimators. 
First we consider Case (B). 



for r^;, = (2r)~^S(X<*,6'*)[r/®2]. Here W is a two-sided standard Wiener 
process independent of Xt* . 

Theorem 2. Suppose that the limit t] = linin-^oo "^n^ {(^l — (^o) exists. Suppose 
that[H]2, [C] and[B] are fulfilled in Case (B) . Thenn'df^{tn—t*)^'^a,Tgmm^^^M.{v) 
as n oo. 

We will prove Theorem [2] and assume for a while that T = 1 to simplify 
the notation. Introduce a new parameter v as t = tl:= t* + v{n^l)-\ Let 



= {M„(4; ^0, ^i) - M„(4; 91)} + {A„(4; ^o, ^i) - An(4; ^t)} 

+{pn{tl9o,9i)-pn{tle;,ei)}. 



Let 
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Lemma 3. For every L > 0, 

sup \D^{v)\^PO 

ve[-L,L] 

as n oo. 

Proof. We assume that v > 0. We have 

1 

^ndeMnii; 6* + n(^o - ^o), ^1 + " ^i*)) du [^;,\9o - 9^ - e{)]. 



For A; = 0, 1 and j = 1, 2, 



sup imM^t- 60,61] 

tG[t*,i*+L{n,?2)-i] 



< [Id^Mr^if + L{n^lr'; 60, 6,)\^] +0(1) 

< 8M-'snpE[\dl^g,{6,)\']+0{l) 

i>l 

Then Sobolev's inequahty imphes 

sup \deMn{t;6o,6i)\ = Op{l). 

te[t* ,t*+L(ni}'^)-''-], 

As a result, 

sup |M„(4;^o,^i)-M„(4;^*,e*)| 

as ^ 00. 

Set r„ = 1^0 — ^ol + 1^1 ~ ^1 1- Simple calculus yields 

|{Tr y - logdet(/d + y)} - {Tr X - logdet(/d + x)}| < csly - x\{\x\ + \y - x\) 

for d X (i-symmetric matrices x and ?/ whenever |a;|, |y| < Cg, where Cg and 
C3 are some positive constants independent of x, y. Indeed, the formula 
/ exp(-2~^(/(i + ex)[z'^'^])dz = (27r)'^/^ det(/rf + ex)"^/^ is convenient for ex- 
plicite computation. 
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Applying this inequality to y = Si^i{9o) ^/^5'i„i(6'i)S'i_i(^o) - h and 
X = Si-i{9Q)~^^'^Si-i{6l)Si-i{6Q)~^^'^ — Id, we see that there exists a constant 
C4 such that for large n, on Bk H {16*^ — 6*\ < (A; = 0, 1)}, 



[nt] 



i=[nt*]+l 

Thereofore, for any e > 0, if we take sufficiently large K, then 



hm„^oo-P 



This implies 



sup \A^{t;9o,9,)-An{t;e*,ei)\>e 
te[f,t*+LM2)-i] 



sup |A„(4;^o,^i)-^n(4;^S,^ni ^"0 



< e. 



as n ^ 00. The convergence 

sup |p„(4;^o,^i)-Pn(4;^0,^l*)l^"0 
t)G[0,L] 

can be shown in the same way as (jS]). 

A similar proof of the uniform convergence on [—L, 0] is possible. After 



all, we obtained the desired result. 



□ 



Remark 4. When 61 {k = 0, 1) are known, we do not need Lemma [3l 

Thus we can focus only on "^nitl; 9^,91) - '^nit*;9Q,9l). For simphcity, 
we write ^'Ji(t) for "^nit; 9^,91). By assumption, there exists a limit t] = 
lim^^oo 'f^n^(^i ~ ^o)- ® denotes the D-space on an interval of t. Let 

and 

iF{v) = -2 (rir w{v) - ^rr,M^ 

for r,,, = l{.>o}(2T)-iH(X,-„r)[<]_ 
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Lemma 4. Let r] = lim„_oo — 9q). Suppose that [H]2, [C] and [B] 
are fulfilled m Case (B). Then H„ ^'^-^^t) D([-L,L]) as n ^ oo for 

every L > 0. 

Proof. We will only consider positive v since the argument is essentially 
the same for negative v. Let T = 1 as before. It follows from Lemma [T] that 



where 



M4t* + vin^l)~';e;,ei), 
A^{t* + v{n^l)-'-e*,ei), 

Pn{f + v{n^lr';9*,9l). 



The evaluation of these terms will be done in the following. As repeated 
previously, we may proceed discussion on the event hereafter. First 



i=lnt*]+l 
fti 



h~'({ r a{xiei)dWtf'-EZ,[ s{x;,ei)dt] 



+Op(l) 



(6) 



where Un{v 



Op(l) means that sup^g[o,L] l^nl"^')! 0> used the 

hypothesis ni)'^ — > oo and the fact that |5'i_i(6'o)~^ — Si-i{9iy^\ < C-dn with 
the localization. To obtain Op(l), L^-estimate helps. It follows from [Hjj 
(i)(c) and (ii) that 



s{xi,ei)dt-El, 



s{xi, ei)dt 



\s{xix,)-s{xi_^,ei)]dt 



ti-l 



[S{Xl,9l)-S{Xl_^,9l)]dt 



- ti-l 

< Cw[o,T]{n-\xy 
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Moreover, with the Burkholder-Davis-Gundy inequahty, the first terms on 
the right-hand side of (El) equals M^{v)+dp{l) with M^{v) = ESSi''^ ^n,i, 



where 



l{r>o}Tr 



■{h~\A,Wf^-Ir) 



(7) 



Now we introduce the backward approximation 



l{r>o}Tr 



ih-\A,Wf'-Ir) 



to ^n,i for en = 2Ln-H-^. After all, by 



^^=[nt*]li^ in,i, we have 
(8) 



Since 



ii.>o}2^;\| V(x;_,„, ei) {s{xj,_,^, e*)-' - s{xj,_,^^,ei)-')a{xj,^,^, ei] 
i{^^o}2^-'vEix^,__^^,e*)m - eX'] + 0,(1) 



the central limit theorem ensures the convergence —2T^^r W in 

D([0, L]). Indeed, the joint convergence of ^;iV(X;_,^, (5(X,-,_,,^, 

S{Xl,^^^,ei)-')aiXl,_^^^,9l) and the process ^„ El=[ntSi'^ {h-\AWr' - 
Ij.) implies the joint convergence of {X^,, M^). In the same fashion, we can 

show -2T\rW{- + L) in D([-L, 0]) if is defined in a natural way 

over negative where {W' {u))ue[o,L] is a standard Wiener process indepen- 
dent of (>V(w))„g[o,L] and J^. Since E{Xl, , e*)[ri®^] is independent of W, we 
can replace the stochastic integral with respect to W in the representation 
of the limit distribution of by the one with respect to the negative-time 
part of the two sided Wiener process W reversible in time. Easy calcula- 
tions yield sup„g[_i ^L] \ A^{v) - r^,^w| -^p and sup^^f.^^^j for 
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extended and p,^ to [—L,L], which completes the proof. □ 

Proof of Theorem We have supposed that T = 1 to state the lemmas, 
and we start with this case. Write v = argmin^g]gE[(t;). For e > 0, take large 
K so that P[t = T] > 1 — e. It follows from Lemma H] that for every a; G M, 

]hRn^^P[n^l{L -t*)<x]-e 

< h^n^ooPl inf Ml{v)< inf MHv)] + snp P[n^l{L - t*) ^ [-L, L]} 

v€l~L,x] ve[x,L] n 

= P[ inf ir{v)< inf IT (v)] + snp P[n^l{in - t*) ^[~L,L]} 

< € + P[v <x]+ P[v ^ [-L, L]] + supP[n<(4 - t*) ^ [-L, L]] 

n 

As L ^ cxD, the last two terms of the right-hand side of the above inequality 
tend to thanks to Theorem [1] (b). So we have obtained 

W^n^^P[n^l{tn -t*)<x] < P[V<X\. 

The estimate of P[nt?^(t.„ — t*) < x] from below can be done in a similar 
manner, which concludes the proof in case T = 1. 

For general T, we introduece a stochastic basis B = {Q, JF, F, P) with 
F = {J^Tu)uelo,i], and the processes = bxu, = ^Tu and = Ym, 
u G [0, 1], to scale the time as t = Tu. Those stochastic processes satisfy the 
stochastic integral equation 

Jo Jo 

where a{x, 9) = y/Ta{x, 6) and W is an r-dimensional F- Wiener process. 
The sampling times {iT/n)2^Q now chage to (z/'n.)"=o setting after 

scaling time. For the change point estimator u„ for u* = T^^t*, we know 

n^l{un-u*) ^'^^ argmin^gigH({}), (9) 

where m{v) = -2{f^W{v) - 2"if = 2-^S(X„., r)[r/®2] and W is a 

two-sided Wiener process independent of cr(X„. ,6^*) = \^a{Xt*,6*). Since 

T argmin^gK m{v) = argmin^g^ ^ (^) 

argmin^gK Il{v) 
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thanks to W(-) T^/'^W{-/T). Thus ^ gives the desired convergence of 
tn since t„ = Tun- □ 

Let us investigate the hmit distribution of the estimator in Case (A). By 
nature of the samphng scheme, only the set Qn = {kT/n; /c G Z} has essential 
meaning for the optimization with respect to the parameter t. Without loss 
of generality, we modify t„ so that it takes values in and set kn = ntn/T. 
Let 

'a{Xt*,oi)is{x,.,e*r' - six,,,ei)-')a{x,.,ei)cr' 

-logdet (^SiXt^,9*)-'S{Xt.,ei) 

where Q are independent r-dimensional standard normal variables indepen- 
dent of Xt* . 

Theorem 3. Suppose that [H]i, [C] and [A] are fulfilled in Case (A). Then 
kfi — [^] argmiUjjg^ K{v) as n ^ oo. 

Proof. We change the definition of t^ and newly set tl = [^]^ + 
Lemma [3] is still valid by essentially the same proof and hence we may only 
consider 6*, 91) - ^*). Writing for 91), we will 

investigate the behavior of the random field 

K„(^;) = ^:iti) - ^:{t*) 

on w e Z. For a while, we consider nonnegative v. The argument is similar 
for negative v. According to Lemma [H we have the decomposition 

Kn{v) = Mn{v) + An{v) + ^„(w), 

where M^iv) = M^{tl; 9*, 9*,), A^iv) = A„(4; 9*,) and Qn{v) = p„(4; ^o, ^D- 
Now, Wln{v) admits a similar expansion as before: 

[nt*/T]+v 

Mn{v) = J2 U^ + Op{l) 

i=[nt*/T]+l 

with ^n.i given by ([7]). Moreover, for e„ = n~^/^ this time, we consider the 
backward approximation of ^n,i, that is. 
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Here v G [0, L] flZ, however this approximation is available when we consider 
V G [—L^ 0]. let Lq be the maximum integer in [0, L\. By continuity of a and 
because W is an F- Wiener process, we have 

where Q are independent r-dimensional standard normal variables indepen- 
dent of XI*] we use the same symbol Q as in the statement. Consequently, 

where 

[nf/T]+v 
i=[nf/T]+l 

and C.oo,i is given by 



l{^>o}Tr 



For A„, we have An{v) Aoo(f) with 

Aoo(^) = i{.>o} Yl \T^[Six;„e*,r'six;„9i)-i, 

i=lnt*/T]+l ^ ^ 



-logdet i^s{x^„e;)-'s{xi,,9i] 

On the other hand, gn{v) tends to uniformly in v. Therefore, 

where ]K^(t>) = Moo('y) + Aoo(f )• Removing r by letting K — > oo, and using 
Theorem [1], we obtain the limit distribution of t^. □ 
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5 Initial estimator for 6k 



In this section, we will breifly discuss the consturction of the initial esti- 
mators. There are two situations according to the prior knowledge of the 
parameter space T of the change point. The first one is the case where 
T = [to,ti] C (0,1) for given numbers to and ti. In the second case, we 
do not assume a prior information of to and ti, instead the precision of the 
initial estimator will be lost. Let 

[nt/T] n 

*n(i;^o)= 5^ G,(^o) and $^(t;^i)= ^^(^i)" 

«=1 i=[nt/T]+l 

Suppose that to and ti are known. Let ^o and 6i satisfy 

^ni^k-Jk) =min<I)^(tfc;^fc) 

Ok 

for k = 0,1. To validate asymptotic properties of the estimators, it is suffi- 
cient that these relations are satisfied asymptotically. Under suitable regu- 
larity conditions as well as the identifiability conditions that 

to rT 

Q{Xt,9*,9) dt > a.s. and / Q{Xt,9*,9) dt > a.s. (10) 

for every 9 ^ 9*, it is possible to show that 9k — 91 = Op{n~^^'^), therefore Con- 
dition [C] is satisfied in both cases (A) and (B). Based on 9k, the estimator t„ 
are defined. According to the previous sections, t„ possesses nt^^-consistency 
and the asymptotic distribution in each case is already known. 

We can also construct the second stage estimators. Let 6„ be a sequence 
of positive numbers such that 6„ = (wd^)'^, where S E (2, oo) is a constant 
satisfying wdf^ — oo as n — > oo. Construct 9k so that 

$^(4 + (-l)'+'6n; 9k) = min <D^(4 + (-1)'+'^^; 9k) 

for k = 0, 1. The new estimators 9k are expected to improve 9k since they 
utilize up to the data near t*. Further, it is possible to construct a new 
change-point estimator with those estimators. Based on 9k, we define t„ for 
t* as 

tn = argmin^g[oy]$„(t;^o,^i)- 
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Since it is usually easy to verify Condition [C] for ^fc, we will then obtain the 
same asympotic results for in as 

Next, let us consider the second situation. The knowledge of is not 
available and it means that any data set sampled over a fixed time interval 
[0, a] is usuless for estimating Oq since t* may be less than a and then the data 
over {t*, a] causes bias in general. A similar notice is also for the estimation 
of 6i. This consideration suggests the use of estimators dk based on the 
data over time interval [0, a„] for k = and the one over [T — a„, T] for 
k = 1, respectively, for some sequence a„ tending to zero. We assume that 
there exitst a constant (3 G (0,1/2) such that a„ > l/{n^l/'^) and that 
\0k — ^fcl = Op((na„)~^) for k = 0,1. When lim^^oo^^n > 0, we also assume 
nan oo. In particular, the first condition implies ni?^ oo. The second 
condition is natural because the number of data is proportional to na„. To 
obtain 6k, we may need the identifiability condition that a{6,x) = a{6',x) 
implies 6 = 9'; it is a strong condition like monotonicity of (j{9, x) in 9. Under 
the assumptions, [C] holds and after that it is possible to construct t„, 9k 
and in in turn as mentioned above. The asymptotic properties of in are the 
same as tn because ^^'s satisfy Condition [C]. It is expected that the new 
estimator in posseses equal or better precision than t„ as numerical studies 
in Section [H] suggest. 



6 Numerical studies 

In this section we run some simulation experiments to asses the quality of 
the estimator of the change point and of the volatilities, under two different 
models. We first consider the following diffusion model without drift 

f Xo + /*(! + XyUWs for t G [0, t*) 
Xt = \ (11) 

[ Xt* + £ (1 + X^fUWs for te[t*,T]. 

where t* is the true change point assumed to be t* = 0.6. The true value of 
the parameters are 9q = 0.2 and 9^ = 9q + n~^, with 7 = i, is the sample 
size and T = nh = 1. The initial value is Xq assumed to be constant, in 
particular we take Xq = 5. The sequences an = bn = with 5 = 3 so that 
they satisfy the properties required in Section [51 The lirst stage estimator 
of 9q (resp. 91) is obtained using the first nan observations from the left 
(resp. nan from the right). We denote the first stage estimators with 9i, 
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i = 0,1. Once the first stage estimators of and 6^ are available, the first 
stage estimator of t*, i.e. t„ is obtained via 

*n(4;^o,^i) = min <l>„(^;^o,^l)• 
^e[o,T] 

Then, with the first stage estimator of t* in hands, we calculate the second 
stage estimator of 9* using observations in the interval [0, t„ — &„] for 9q and 
observations in the interval [t„ + bn,T] for O^. We denote the second stage 
estimators of 9* by 9i. Finally, the second stage estimator of t*, i.e. in, is 
obtained as 

*n(tn;^o,^i) = mm <!>n{t;9o,9i). 

ie[0,T] 

For comparison, we also report the value of the estimator t„ obtained 
plugging the true parameter values in the contrast function, i.e. when the 
volatilities are supposed to be known 

<^n{in,9;,9l) = min <^n{t;9*,9l), 

iG[0,T] 

and this can be considered as a benchmark. For the Monte Carlo setup, 
we consider different sample sizes n = 1000, 2000, 5000 and for each sample 
size n, we run M = 10000 Monte Carlo replications. Under this choice of 
n the value of 91 = 0.3778, 0.3495, and 0.3189 respectively. The values of 
the sequences a„ and 6„ are reported in Table [TJ Observations are supposed 
to be sampled at sample rate h = 1/n. Table [T] also reports Monte Carlo 
estimates (i.e. average over the M replications) of the volatility parameters 
6*0 and 9i and the change point t*. In parenthesis are the standard deviations 
of the Monte Carlo estimates. In the second experiment we consider a Cox- 
Ingersoll-Ross (1985) model 

r Xo + /J ,/9*X:sdW, for t e [0, t*) 
Xt={ , (12) 

I Xt, + ,/9lX:dWs for t e [t*, T]. 

with change point t* = 0.7 and all remaining experimental conditions are the 
same as in previous experiment. The results are reported in Table [2l The 
difference in the two experiments is only in the regularity of the diffusion 
coefficient term. Comparing the two simulation results, it is possible to see 
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that the second stage estimators in the second experiment performs shghtly 
better in term of the standard deviation. 

We also consider the behaviour of the asymptotic distribution of the 
change point estimator for second stage estimator in the first model, for 
sample size n = 5000. In particular, due to mixed-normal limit, we studied 
the distribution of the studentized limiting distribution of n^^(t„ — t*) under 
the true model, i.e. 

z = neliL-t*)T{Xt^,eo), 

with t{Xf>,9o) = (log(l + X2))2. Then Z converges to W{v) - l\v\ with 
density 



and distribution function 

Fix) 




X > 
x<0 



with the distribution function of the gaussian random variable, and 



9{x) = 1 + y l^e-f - i(x + 5)$ + ^e^'$ (-^^S 

(see e.g. Csorgo and Horvath, 1997). In Figure [1] we report the graphical 
representation of the histogram and empirical distribution function of Z (over 
10000 Monte Carlo replications) against their theoretical counterparts which 
looks quite reasonable. 
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histogram vs limit distribution (second stage) 
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EIDF vs limit (second stage) 
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Figure 1: Histogram versus theoretical density function (up) and empirical 
distribution function versus theoretical distribution function (bottom) for 
the second stage change point estimator. Results of 10000 Monte Carlo 
replications and sample size n = 5000 for the first model. 
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n 


an 




eo 


01 


tn 


^0 


ei 


tn 


5000 


0.1189 


0.601 
(0.005) 


0.200 
(0.009) 


0.319 
(0.014) 


0.601 
(0.011) 


0.200 
(0.005) 


0.319 
(0.013) 


0.601 
(0.012) 


2000 


0.1495 


0.601 
(0.008) 


0.200 
(0.013) 


0.349 
(0.020) 


0.601 
(0.014) 


0.200 
(0.008) 


0.349 
(0.017) 


0.601 
(0.015) 


1000 


0.1778 


0.601 
(0.011) 


0.199 
(0.017) 


0.377 
(0.025) 


0.601 
(0.019) 


0.200 
(0.011) 


0.377 
(0.026) 


0.602 
(0.018) 



Table 1: Monte Carlo estimates for model (|11|) over 10000 replications. True 
values: 9* = 0.2, 91 = 0.378, 0.350, and 0.319 for different sample sizes n = 1000, 
2000 and 5000. True change point t* = 0.6. 



n 


an 


tn 


00 


9i 


in 


^0 


9i 


tn 


5000 


0.1189 


0.701 
(0.010) 


0.200 
(0.012) 


0.319 
(0.018) 


0.701 
(0.011) 


0.200 
(0.018) 


0.319 
(0.012) 


0.701 
(0.010) 


2000 


0.1495 


0.702 
(0.016) 


0.200 
(0.016) 


0.350 
(0.029) 


0.701 

(0.024) 


0.200 
(0.009) 


0.350 
(0.030) 


0.701 
(0.021) 


1000 


0.1778 


0.703 
(0.025) 


0.200 
(0.021) 


0.378 
(0.040) 


0.701 
(0.038) 


0.200 
(0.012) 


0.377 
(0.056) 


0.701 
(0.040) 



Table 2: Monte Carlo estimates for model ([T2]) over 10000 replications. True 
values: 9^ = 0.2, 0* = 0.378, 0.350, and 0.319 for different sample sizes n = 1000, 
2000 and 5000. True change point t* = 0.7. 
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